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Abstract 

Simulations of the two-dimensional self-avoiding walk (SAW) are performed in a half- 
plane and a cut-plane (the complex plane with the positive real axis removed) using the 
pivot algorithm. We test the conjecture of Lawler, Schramm and Werner that the scaling 
limit of the two-dimensional SAW is given by Schramm's stochastic Loewner evolution 
(SLE). The agreement is found to be excellent. The simulations also test the confor- 
mal invariance of the SAW since conformal invariance implies that if we map infinite 
length walks in the cut-plane into the half plane using the conformal map z — > y^, then 
the resulting walks will have the same distribution as the SAW in the half plane. The 
simulations show excellent agreement between the distributions. 
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1 Introduction 



Lawler, Schramm and Werner |2j have conjectured that the scahng hmit of the two-dimensional 
self-avoiding walk (SAW) is given by Schramm's stochastic Loewner evolution (SLE). SLE 
is a two dimensional conformally invariant random process which depends on a parameter k, 
and so is denoted SLE^. Chordal SLE refers to the version of SLE in which the random curve 
or set connects two points on the boundary of a simply connected domain. It is usually defined 
first for the case where the domain is the half-plane and the two boundary points are and oo. 
Its definition is them extended to other simply connected domains D and boundary points using 
a conformal map from the half-plane to D which maps the two boundary points appropriately. 
If K < 4, chordal SLE gives a probability measure on simple curves, i.e., curves that do not 
intersect themselves pilj. The conjecture of Lawler, Schramm and Werner is that for any simply 
connected domain D and points z and w on its boundary, SLEg/s is the scaling limit of SAW's 
that go from z to w and stay inside D. 

For K = 8/3, Lawler, Schramm and Werner |HI have a theorem that makes it possible to 
explicitly compute the distributions of many random variables associated with the SLE random 
curve. For the scaling limit of the SAW, these random variables can be studied by simulation. 
Thus one can numerically test their conjecture that the scaling limit of the SAW is SLEg/a 
by comparing the distributions from simulations of the SAW with the exact distributions for 
SLEg/3. This test was carried out for two such random variables for the SAW in the upper 
half-plane in 0, and excellent agreement was found. In this paper we consider more random 
variables for which the exact distribution can be computed for SLEg/3. We compare their exact 
distributions with the numerical distributions of the same random variables for the SAW in the 
half-plane. We also simulate the SAW in the cut-plane consisting of the complex plane minus 
the non-negative real axis. The map z — ^/z takes the cut-plane onto the half-plane, and 
by composing the random variables for the half-plane with this map we obtain corresponding 
random variables for the cut-plane. We compare their distributions for the SAW from our 
simulations for the cut-plane with the exact distributions for SLEg/3. We also consider the 
probability that the walk passes to the right of a given point in the half-plane (or the cut-plane) 
and compare this probability for the SAW simulations with an exact formula of Schramm ^2] 
for SLE. 

Note that for both of the domains we consider, the terminal point of the walk is at infinity. 
This case is particularly well suited to simulations, since it is expected that we can construct 
the scaling limit by considering all SAW walks with a fixed length which start at the origin, 
taking the limit N ^ 00 and then taking the limit that the lattice spacing goes to zero. For a 
domain D and two finite points z and w on its boundary, the scaling limit should be constructed 
as follows. We introduce a lattice and consider all self-avoiding walks which start at z and end 
at w. The probability of such a walk is taken to be proportional to (3~^ where N is the number 
of steps in the walk, and j3 is the constant such that the number of SAW's in the plane starting 
at the origin grows with the number of steps, N, as The measure is normalized so that it 
is a probability measure. We then take the limit of this measure as the lattice spacing goes to 
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zero. The construction of the scahng hmit in the case of SAW's with infinite terminal point is 
rather different from the case of a finite terminal point, so it would be interesting to test the 
conjecture that the scaling limit is given by SLEg/s in the case of a finite terminal point. 

In addition to describing the scaling limit of the SAW, SLE is conjectured to describe the 
scaling limit of a large number of other two dimensional models. Many of these conjectures have 
been proved recently. Schramm showed that if the loop-erased random walk has a conformally 
invariant scaling limit, then that limit must be SLE2 [TT]. He also conjectured that the scaling 
limit of percolation should be related to SLEg, and the scaling limit of uniform spanning trees 
(UST) is described by SLE2 and SLEg. The conjectures for the loop-erased random walk and 
the UST have been proved by Lawler, Schramm and Werner j^]. Smirnov has proved the 
conformal invariance conjecture for critical percolation on the triangular lattice and that SLEg 
describes the limit f^]- Lawler, Schramm and Werner used SLEg to rigorously determine the 
intersection exponents for Brownian motion and proved a conjecture of Mandelbrot that the 
outer boundary of a Brownian path has Hausdorff dimension 4/3 jSllllEl- The random cluster 
representation of the Potts model for < g < 4 was conjectured by Rohde and Schramm to be 
related to the SLE process as well J^. 



2 SLE predictions 

The random variables we consider are defined for curves in the upper half-plane as follows. 
Note that these random variables are defined both for the SAW and for SLE. We use 7 to 
denote the random curve in both cases. Consider a horizontal line at a height of c above the 
horizontal axis. The curve 7 will intersect it, possibly more than once, and we look for the 
left-most intersection. The random variable is the x-coordinate of this intersection, divided 
by c. So 

Xe = - min{a; : x + ic G 7} (1) 

c 

We can also consider the first intersection of the curve with the horizontal line. ("First" means 
the first intersection as we traverse the curve starting at the origin.) We let Xf be the x- 
coordinate of this intersection, divided by c. (The subscripts e and / are for "extreme" and 
"first," respectively.) The next random variable is defined using a vertical line at a distance c 
to the right of the origin. The curve will intersect it, and we look for the lowest intersection. 
The random variable Ye is the y-coordinate of this intersection, divided by c. So 

Ye = - mm{y : c + iy E 'j} (2) 

c 

The random variable Yy is the y-coordinate of the first intersection, divided by c. For the final 
random variable, consider a semi-circle of radius c centered at the the point cd on the real 
axis where \d\ < 1. So the origin where the random curve starts is inside the semicircle. The 
intersections of the random curve with the semicircle are of the form c{d + e*^) and we look for 
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the intersection with the smallest 6. The random variable Ge is this smallest angle, normalized 
so that it ranges between and 1. So 

Ge = - min{^ : c{d + e'^) G 7} (3) 

TT 

The random variable G/ is the angle of the first intersection, again normalized so that it ranges 
from to 1. If the probability measure is invariant under dilations, then the distributions of 
all of these random variables are independent of c. This is true for SLE and is expected to be 
true for the scaling limit of the SAW. 

The distributions of Xe, Y^-, Ge are all easily computed using the following theorem of Lawler, 
Schramm and Werner. Let EI be the upper half-plane. Let A be a compact subset of the closure 
of EI such that EI \ A is simply connected and is not in A. Let ^ be the conformal map from 
EI \ y4 onto EI which fixes and 00 and has $^(00) = 1. We continue to denote the random 
curve generated by SLE, the SLE "trace," by 7. 

Theorem 1. (Lawler, Schramm, Werner 181.) For k, = 8/3, chordal SLE in the half plane has 

P{^nA = (ll) = (!?'^{Of (4) 
Our next step is to use this theorem to compute the distributions of Xe, Ye and Gg. 



2.1 Hitting the horizontal hne 

It is convenient to take c = vr to compute the distribution of X^. Let Lt be the horizontal ray 
which starts at t + ni and goes to the left. Let (z) be the conformal map which maps H \ 
onto EI and satisfies the conditions in the theorem. Note that Xg < t/vr if and only if 7 hits Lt. 
So by the theorem 

PiXe<t/n) = 1-^1(0) (5) 

The map w{z) = z + ln{z) + 1 + t maps EI onto EI \ L^. We need the inverse of this map but it 
cannot be explicitly found. The inverse should be normalized so that it fixes and 00 and has 
derivative 1 at 00. The above map does not fix 0, but meets the other two conditions. Fixing 
is not necessary since we can achieve this condition by just adding a constant to the inverse 
map, and this which will not change its derivative. So we have 

where zq is the image of under the inverse map, i.e., = + ^^i^o) + ^ + t. 
Define 

g{x) = X + ln(x) (7) 
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Figure 1: The distribution, P{Xe < t), of for the half-plane. The solid line is the distribution 
for SLEg/s, and the open circles are the results of the simulation of the SAW. 



This is an increasing function which maps (0, oo) onto the real line, so it has an inverse that 
maps the real line to (0, oo). Note that zq = g^^{—t — 1). We have 

dw , , 1 
So (0) and a trivial change of variables gives 

Although cannot be explicitly computed, it can be trivially computed numerically. The 
graph of the above distribution is the solid line in figure ^ The open circles in the figure are 
the results of the simulation for the SAW. 

We can find the asymptotic behavior of the distribution in Q as t goes to ±oo. For large 
positive t, g{t) = t + ln(t) ~ t. So as t ^ — oo, g~^{—7Tt — 1) ^ —nt, and so 

/ -nt 5 

as i^-oo (10) 
As t 0, g{t) ^ ln(t). So as t ^ oo, g'^i-nt - 1) ^ e^^^-\ So 

/ „-7rt-l \ 5/8 

P(X,<t)^l-(-^^-^ ^ 1 - e-^(-*+^)^ a. t->oo (11) 
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Figure 2: The distribution of for the half-plane. The solid line is SLEg/s, and the open 
circles are the SAW. 



As t — oo, the probability goes to zero slowly, but as t oo, the probability goes to one 
exponentially fast. This is reasonable since when Xe is very negative it only means there is at 
least one intersection with the horizontal line far to the left of the origin, but when Xg is very 
positive it means that all intersections with the horizontal line are far to the right of the origin. 

2.2 Hitting the vertical line 

The distribution of was studied in We take At to be the line segment from 1 to 1 + it. 
The conformal map that maps HI \ onto EI with the required normalizations is 

^^^{z) = t,/-{z-iy-t^ (12) 

where the square root has its branch cut along the negative real axis. Thus the distribution of 
Ye is 

P{Ye <t) = P(7[0, oo) n A 7^ 0) = 1 - '^'^,(0)'/' = 1 - (1 + t')~'/'' (13) 
Figure |21 shows this distribution and the results of the simulation for the SAW. 

2.3 Hitting the circle 

It is convenient to translate so the semicircle is centered at the origin. Setting c = 1, this means 
the random curves start at -d. So P{Qe < t) = 1 - $'(-rf)^/^ where $ is a conformal map 
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Figure 3: The distribution of for the half-plane for d = 0, 0.5, 0.9. (d increases from left to 
right.) The solid lines are SLEg/s, and the open circles are the SAW. 



which takes the half-plane minus the arc = {e*^ : < ^ < irt} onto the half-plane with the 
normalizations that the map fixes oo and has derivative 1 at oo. (As in the previous case, we 
ignore the condition that the map fixes the origin since it does not affect the derivative.) 
The conformal map 

z ^ ^ (14) 

sends the upper half-plane (including oo) onto itself, and it sends the upper half of the unit 
circle to the upper half of the imaginary axis. Let 



a — 



sin(7rt) 
1 -|- cos(7rt) 



(15) 



Then the arc At is mapped onto the line segment from to ia. We can then map EI with this 
line segment removed onto EI as we did in the previous section. Composing these two maps 
and multiplying by a factor of vT+~a^ for later convenience, we define 



f^-11 



(^ + 1)^ 



1/2 



(16) 



with the branch cut for the square root being the negative real axis. The map i) sends H \ 
onto H . 
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The map ip does not send oo to itself. For z near oo, 



^(z) = (l + a^)(l 
In particular, 7/^(00) = (1 + a^). Now let 



+ 



1 + a2 - iIj{z) 



For large z, ^Ati^) ~ z, so the derivative at 00 is 1 as required. 
For real x with — 1 < x < 1, the choice of branch cut leads to 



= -VlT 



Define 



2\-l 



(a; + 1)2 
1 + cos(7rt) 



This will prove to be a natural variable to use. We have 



^{x) = j= 



(x + l) 



nl/2 



- - 1 



SO 



1 

s 
2s 



1 - 



Axs 



.1/2 



(x + l) 



Computing the derivative $'(— (i) then yields 

PiQe <t) = l- 



1 + 



Axs 



1/2 



4s2(l + d) 



1 - d + [(1 - d)2 + Adsyl^Y ((1 - + Ads) 



1/2 



(17) 



(19) 



(20) 



{2V. 



[22) 



(23) 



For c? = 0, 0.5 and 0.9, this distribution and the results of the simulation for the SAW are shown 
in figure El 



2.4 Passing right 

In addition to the distributions of the random variables Xg, Y^. and Gg, we also consider the 
following probability. Fix a point in the upper half-plane. One can then ask if the random 
curve passes to the right or left of this point. For SLE this probability only depends on the 
polar angle of the point since SLE is invariant under dilations. This should also be true for 
the scaling limit of the SAW, since it is expected to be invariant under dilations. Schramm 
P2] rigorously derived an explicit formula for this probability for k < 8. For general n it is 
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Figure 4: The probability that the walk passes to the right of a point as function of its polar 
angle for walks in the half-plane. The horizontal axis is the angle divided by vr, so that it ranges 
from to 1. The solid line is the exact result for SLEg/s, and the open circles are the results of 
the simulation of the SAW. 



given by a hypergeometric function, but for k = 8/3, his formula is quite simple. Denoting the 
probability that the curve passes to the right of a point with polar angle 6 by p(6'), he showed 
that for K = 8/3 

p{e) = \{i-cos{e)) (24) 

In our simulations we study this probability by fixing a radius c and computing the probability 
the path passes to the right of ce*^ for a large set of values of 9. The above function and the 
results of the SAW simulation are shown in figure El (Note that the horizontal axis in the figure 
is e/Ti.) 

2.5 The cut-plane 

The cut-plane we consider is the plane with the non-negative real axis removed. Let f{z) = 
with the branch cut along the positive real axis. Then / maps the cut-plane onto the upper 
half-plane. We will continue to denote curves in the upper half-plane by 7, and use 7 to denote 
curves in the cut-plane. Given a curve 7 in the cut-plane, 7 = /07 is a curve in the upper 
half-plane. So we can define the various random variables for the cut-plane by applying their 
definitions in the half-plane to /07. We will put a " on top of random variables defined on 
curves in the cut-plane. For the simulations it is useful to work out these definitions explicitly 
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in terms of the curve 7 in the cut-plane, rather than map each SAW in the cut-plane to the 
half-plane. 

First consider ©e and G/ for d = 0. The map / simply divides the polar angle by 2, so for 
curves 7 in the cut-plane, 

Oe - 7^ min{^ : ce*^ e 7} (25) 

The random variable ©/is the polar angle of the first intersection of 7 with the circle, divided 
by 2tt. If d 7^ 0, the image of the semicircle under 2; — > is not a circle. We have not simulated 

Qe or Qf in this case. 

To find the definition of X^, we first take c = 1. The image of the horizontal line {i + 1 : 
— 00 < i < 00} under 2; — > is a parabola whose axis is the horizontal axis and which opens 
to the right, 

x^t^ -1, y^2t (26) 

In the half-plane, Xg is the smallest t such that i -|- i is on the curve. In the cut-plane, we 
consider all intersections of 7 with the parabola and find the intersection with the smallest 
^/-coordinate. Since t — y/2, X^ is one half of the y-coordinate of this "lowest" intersection. 
Equivalently, 

Xe = mm{t : (t^ - 1, 2t) E 7} (27) 

SLE is invariant under dilations of the cut-plane, and the scaling limit of the SAW in the 
cut-plane is expected to have this invariance as well. So for c 7^ 1 we can take the parabola to 
be 

x = c(i^-l), y^2ct (28) 

and let 

Xe = mm{t : {c{t'^ - 1), 2ct) G 7} (29) 

Xf is the y-coordinate of the first intersection of 7 with the parabola divided by 2c. 

To find the definition of 1^, we consider the image of{l-|-it:0<t<oo} under z ^ z^. It 
is the upper half of a parabola whose axis is the horizontal axis and which opens to the left: 

x^l-t"^, y^2t, t>0 (30) 

In the half-plane, Y^. is the smallest t such that 1 -|- it G 7, so in the cut-plane Yg is one half of 
the y-coordinate of the lowest intersection of 7 and the half parabola. 

Fe^niinO : (l-t^2t) e7,i >0} (31) 

More generally, we can let 

Ye = minji : (c(l - t^), 2ct) e 7, t > 0} (32) 
Yf is the y-coordinate of the first intersection with the parabola divided by 2c. 
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We have defined the random variables in the cut-plane so that if the probabihty measure 
is conformally invariant, then they will have the same distribution as their counterparts in the 
half-plane. Rather than compare the distributions of the random variables Xe, Ye and 0e with 
those of Xe, Ye and 6e, we will compare all these distributions with the SLEg/s predictions, eqs. 
(jni), (fT!?jl and ((2S1)- This tests both the conjecture that the scaling limit of the SAW is SLEg/s 
and the conformal invariance of the SAW. For the random variables Xf,Yf and 6/, we do not 
know their distributions for SLEg/s. So we will directly compare the distributions of Xf,Yf 
and 0/ with those of Xf, Yf and O/. This tests the conformal invariance of the SAW. 

3 The simulations 

In all of our simulations the walks had one million steps. For the half-plane we ran the pivot 
algorithm for 10 billion iterations of the Markov chain. For the cut-plane we ran for 11.4 
billion iterations. The simulation of 0e for d 7^ in the half-plane was done separately and 
consisted of 6.8 billion iterations. For walks with a million steps only about 5% of the proposed 
pivots are accepted. Of course, accepted pivots do not produce independent walks and for the 
random variables considered here most accepted pivots do not even change the value of the 
random variables. So the number of effectively independent samples is considerably less than 
the number of accepted pivots. Each of the simulations requires about a month on a 1.5 GHz 
PC. The exact speed of the simulation depends on the choice of the half-plane vs. cut-plane 
and how many random variables are simulated. 

A walk with N steps is typically of size A^^^^, so to study the various random variables we 
take c = IN^/^^ where / is fairly small. Note that if we rescaled to make c equal to 1, the lattice 
spacing would be (/A^^/^)~^. We will refer to this quantity as the "effective lattice spacing." 
Note that I is the ratio of the scale used to define the random variable to the scale of the walk. 
So we must take I small to make the effect of the finite length of our walks negligible. But 
as / gets smaller, the effective lattice spacing gets larger. There is a second effect as / gets 
smaller. For smaller /, the fraction of the pivots that change the values of the random variables 
is smaller. So the statistical errors get larger as / gets smaller. We do not know a priori what 
value of / will be optimal, so we compute the distributions of each random variable for four 
different values of / in our simulations. The particular values of / that we use are determined by 
some experimentation with much shorter simulation runs. We do not use the same four values 
of I for the different random variables. 

In figures Q to 01 we show the distributions of Xe, Ye and Gg. (Throughout this paper 
we work with the cumulative distributions of our random variables rather than their densities 
since any simulation computes cumulative distributions. Computing densities requires taking 
numerical derivatives of the cumulative distributions, and so the densities would have larger 
statistical errors.) The solid curves are the exact distributions for SLEg/3. The circles are the 
results of the simulation of the SAW. Figure El studies the probability that the walk passes to 
the right of a point in the upper half-plane as a function of the polar angle of the point. The 
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Figure 5: Half-plane: The distribution of for the SAW minus the distribution of Xg for 
SLEg/a- The top curve, with the larger error bars drawn with solid lines, has / = 0.01, and the 
bottom curve has I = 0.05. 



solid curve is Schramm's exact result for SLEg/a, and the circles are the results of the SAW 
simulation. In all of figures Q to EJ one cannot see any difference between the SAW simulations 
and the exact curves for SLEg/s. In figures El to 01 we plot the same four quantities, except that 
now we plot the result of the SAW simulation minus the SLEg/3 functions. The first thing that 
should be observed in these figures is the scale of the vertical axis. It is quite small. In all but 
one of these figures the total vertical range shown is 0.007 or 0.7%. In figure (Hlit is 0.008. 

In figures to (HI several values of / are shown. The nonzero effective lattice spacing means 
that we are simulating discrete random variables. So their distributions will be discontinuous. 
After subtracting off the continuous SLE distribution, the jumps will appear in the difference as 
rapid oscillations. As / increases, the effective lattice spacing decreases, and so the oscillations 
are usually "faster" but smaller in amplitude. Also, for larger / a larger fraction of the pivots 
change the values of the random variables, and so a larger I typically produces smaller statistical 
errors. Both of these effects can be seen in all four of the plots. 
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Figure 6: Half-plane: The distribution of for the SAW minus the distribution of Yg for 
SLEg/a. The top curve, with the larger error bars drawn with solid lines, has I — 0.002, and 
the bottom curve has I = 0.005. 
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Figure 7: Half-plane: For d = 0, the distribution of 0e for the SAW minus the distribution of 
©e for SLEg/s. The three curves shown are for / = 0.2,0.1,0.05 (from top to bottom). As / 
decreases the finite length effects decrease, but the error bars and lattice effects grow larger. 
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Figure 9: Half-plane: The probability that the SAW passes to the right of a point as function 
of the polar angle of the point. The corresponding function for SLEg/a has been subtracted off. 
Going from top to bottom on the left half of the figure, the curves are I — 0.2, 0.1, 0.05. 
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As / becomes larger, the effect of the finite length of the walk will begin to be seen. This is 
well illustrated by figure [71 which shows the distribution of 6e for d = 0. For the largest value 
of I shown, I = 0.2, the effect of the finite length of the walk is clear - the curve differs from zero 
by several times the size of the statistical errors. This curve is the smoothest of the three curves 
and has the smallest statistical errors. For Z = 0.1 the finite length effect is greatly reduced, but 
is still statistically significant. The / = 0.05 curve seems to be the best of the values of / that 
were simulated. The maximum difference of the SAW and SLEg/s distributions is only about 
0.05%. Our simulations included a fourth value of / which is not shown, I = 0.02. For this value 
the larger effective lattice spacing and larger statistical errors produce a difference curve that 
is rougher and larger than the / = 0.05 curve. The behavior in figure |H1 for the distribution of 
9e for d = 0.9 is quite similar to figure [71 except that the nonzero / effects appear to be larger. 

In figure [S] the finite length effect is clearly seen in the / = 0.05 curve; for large negative 
values of t the deviation of this curve from zero is caused by the walk being too short. In figure 
[B] there are no obvious finite length effects; the deviation of the curve from zero appears to be 
primarily caused by the nonzero effective lattice spacing. The deviation is of the same order 
as the error bars and the oscillations. In figure [HI the finite length effects and nonzero effective 
lattice spacing effects are similar to those seen in figure [7| Note that the I = 0.2 and Z = 0.1 
curves are significantly different from zero at the right, corresponding to a polar angle of vr. 
This effect is a result of the nonzero probability that the walk does not reach the semi-circle or 
that it crosses it, but ends inside the semi-circle. In both of these cases it is unclear whether 
the walk will pass to the right or left of the points on the semicircle. The algorithm must make 
some arbitrary choices in these cases. 

Figures [TU] through [THI show the same quantities as figures [Sj to [7| and [HI but for the cut-plane. 
For the random variable Ge (figure [T^ and the probability of passing right of a point (figure 
IT^ . the agreement is again excellent. In both of these figures the vertical scale is 0.007, the 
same as in the corresponding figures for the half-plane. For the random variables and Y^, 
figures and ^2 the agreement is not quite as good, but the deviations from the SLE results 
are still small. (In these two figures the vertical scale is two to three times larger than in the 
other figures.) For these two random variables it is harder to do accurate simulations for the 
following reason. In the cut-plane, X^. and depend on the intersections of the random curve 
with parabolas. It typically takes a longer length of curve to attain these intersections than for 
the lines involved in the definition of X^ and Yg in the half-plane. So in the cut-plane we must 
use smaller values of I. For X^ in the cut-plane, the curves shown use I = 0.002 and / = 0.005 
as compared to I = 0.01 and / = 0.05 for X^ in the half-plane. Even with these small values 
of /, the finite length effects are still quite visible in figure ^| The deviation of the curves 
from for the most negative values of t is pronounced. This is the part of the distribution 
that is particularly sensitive to the need for very long walks to hit the parabola. Of course, 
small values of / mean a large effective lattice spacing and large statistical errors. For Yg in the 
cut-plane, the values of / shown are 0.0005 and 0.001, as compared to 0.002 and 0.005 for Ye in 
the half-plane. The finite length effects in figure ^2 can be seen in the substantial deviation of 
the curves from for large t, again a reflection of the need for long walks to reach the parabola. 
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Figure 10: Cut-plane: The distribution of for the SAW minus the distribution of for 
SLEg/s. The top curve, with the larger error bars drawn with solid lines, has I — 0.002, and 
the bottom curve has / — 0.005. 




Figure 11: Cut-plane: The distribution of Y^. for the SAW minus the distribution of Y^. for 
SLEg/a. The top curve, with the larger error bars, has / = 0.0005, and the bottom curve has 
I = 0.001. 
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Figure 12: Cut- plane: The distribution of ©e for the SAW minus the distribution of ©g for 
SLEg/a. The curve with the greater deviation from the horizontal axis and the error bars drawn 
with dashed lines has / = 0.05. The other curve has / = 0.02. 
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Figure 13: Cut-plane: The probability that the walk passes to the right of a point as function 
of the polar angle of the point. Going from top to bottom on the left half of the figure, the 
curves have / = 0.1, 0.05, 0.02. 
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The scaling limits for the SAW in the half and cut-planes are conjectured to be related by 
the conformal transformation, but there is no reason that the finite length effects in the two 
cases should be related. Indeed, the simulations show they are quite different. For example, 
compare the curves for the largest values of I in figures 171 and IT^ The curve in figure [3 is always 
positive, looking roughly like the first half of a sine wave, while the curve in figure El is both 
positive and negative. 

Finally, we consider the random variables Xf, Yf and 0j in the half and cut-planes. We 
don't know the exact distributions of these random variable for SLEg/s, but we can still compare 
the distributions we get from the simulations of the SAW in the half-plane with the simulations 
for the cut-plane. Recall that Xf, Yf and B/ (the random variables in the cut-plane) were 
defined so that they will have the same distribution as their counterparts in the half-plane if 
the SAW is conformally invariant. If we simply plot the distributions themselves, they agree 
so well that the difference cannot be seen in the plots. So instead of plotting the distributions, 
we plot the distributions minus various reference functions. These reference functions are quite 
ad hoc. They are chosen to be simple functions that are relatively good approximations to the 
distributions. They are defined as follows. For Xf and Xf we use the function 

F{t) = ^ (tanh(1.16t) + 1) (33) 

For Yf and Yf we use the distribution of Yg for SLEg/s, i.e., 

F(t) = 1- (l + t2)-Vi6 (34) 

For Qf and 0j we use 

F{t) = t - 0.12 sin(27rt) - 0.009 sin(47rt) (35) 

We emphasize that these are not meant to be highly accurate approximations of the distribu- 
tions of Xf, Yf and 0/. One could find better approximations with more complicated functions. 
The only purpose of these functions is to provide a convenient reference with respect to which 
we can plot the distributions for the half and cut-planes and compare them. 

Figures IT^ to IT^ compare the distributions of Xf, Yf and 0j in the half-plane with their 
analogs for the cut-plane. Again, the most important features of these graphs is the small scale 
of the vertical axis. For Xf and G/ the difference between the distributions in the half and 
cut-planes is very small. For Yf the difference is somewhat larger for large values of t, but still 
small. We attribute this greater difference to the larger finite length effects in the cut-plane. It 
can take a walk in the cut-plane a long time to reach the parabola involved in the definition of 

Yf- 
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Figure 14: For the half and cut planes the distribution of Xf for the SAW simulation minus 
the reference function ()33|) is shown. In the half-plane / = 0.05, and in the cut-plane I = 0.02. 
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Figure 15: For the half and cut planes the distribution of Yf for the SAW simulation minus 
the reference function fl34j) is shown. In the half-plane simulation / = 0.005. In the cut-plane 
simulation / = 0.001. Even with these small values of /, the finite length effects produce a 
noticeable difference between the curves for large t. 
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Figure 16: For the half and cut planes the distribution of G/ for the SAW simulation minus 
the reference function fj35p is shown. The half-plane simulation used / = 0.1 and the cut-plane 
simulation used I = 0.05. The half-plane curve has error bars drawn with solid lines, while the 
cut-plane uses dashed error bars. 
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4 Algorithmic considerations 



The pivot algorithm is used for our simulations. (This algorithm is discussed in p.) The 
algorithm picks a site at random along the walk, called the pivot point, and picks a random 
element of the group of symmetries of the lattice about the point. This group element is applied 
to the part of the walk after the pivot point. The result is a new nearest neighbor walk, but 
it need not be self-avoiding or lie in the upper half-plane. The walk is accepted only if both 
of these conditions are meet. Otherwise the proposed walk is rejected and the current walk 
is counted as another state in the Markov chain. The Markov chain trivially satisfies detailed 
balance. In the appendix we show that it is irreducible. 

The speed of the pivot algorithm is typically measured by considering the average time 
needed to produce an accepted pivot. The algorithm may be implemented so that this 
time grows with the number of steps, N, as 0{N^) with g < 1. The exact value of q is not 
known and probably depends on details of the implementation, but simulations indicate the 
implementation in pp has q < 0.57 in two dimensions. (This estimate is based on simulations 
of the walk in the full plane, not the half or cut planes.) 

There are two main steps in the pivot algorithm, and both would seem to require a time 
0{N) per accepted pivot. The first is the test for self intersections to see if the new walk should 
be accepted. The second is actually carrying out the pivot. To test for self-intersections quickly, 
we take advantage of the fact that the walk u only takes nearest neighbor steps. Rather than 
simply checking if uj{i) = uj{j), we compute the distance d = \\uj{i) — u;(j)||i. If d is nonzero 
then we can conclude not just that uj{i) ^ but also that 

u;(0^^(/), ^f \i-i'\ + \j-f\<d (36) 

Thus we can rule out a large number of potential self intersections if d is large. Since it takes 
a time 0{N) to simply write down a walk with steps, the second step of carrying out the 
pivot would seem to require a time that is 0{N) per accepted pivot. To do better, the key 
idea is to not carry out the pivot each time a pivot is accepted. Instead we keep track of which 
pivots have been accepted and only carry them out after a certain number have been accepted. 
Details of this implementation of the pivot algorithm may be found in p]. 

In the usual implementation of the pivot algorithm one chooses the pivot point by giving 
equal probability to all the points on the walk. One can, however, take the probability of 
picking the ith site along the walk to be p{i), where p{i) is a function whose sum is 1. The 
only constraint is that the p{i) must be positive. If one is interested in the distribution of 
the end-point of the walk, then every accepted pivot changes this random variable. For this 
random variable it does not appear that anything could be gained by making p{i) non-uniform. 
However, there is a substantial benefit to using a non-uniform p{i) for the random variables 
in this paper. All of our random variables typically depend only on a short segment of the 
walk near the origin. (The smaller / is, the shorter the segment.) So most accepted pivots do 
not produce any change in the random variable. This suggests that it might be worthwhile to 
choose pivot locations near the start of the walk more often than pivot locations far from the 
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start. For the simulations in this paper we define p{i) as follows 



p{i) 



8c, 
4c, 
2c, 



C- 



if < z < |iV 
if^N <i<lN 
if < i < fiV 
if < i < 

5 — 



(37) 



where c = jqN~^ so that the sum of the p{i) is 1. This is a rather ad hoc choice, but a crude 
test indicates that for a given number of iterations of the algorithm, it typically reduces the 
standard deviation of the random variable by a factor of two. A systematic study of the effect 
of p{i) would be useful. 

For each of the six random variables we consider four different values of /. We also consider 
four values of / for the probability of passing to the right of a given point. Thus there are 
28 different observables to be computed, and some care is necessary to be sure that the time 
required for this part of the simulation does not dominate the simulation. All of these observ- 
ables require finding intersections of the walk with a given curve (a line, parabola or circle). 
Searching through the walk one step at a time for these intersections would be disastrous, since 
it would require a time 0{N). Such a search is easily avoided. At a given site in the walk 
we do not simply check if the next step intersects the curve. Instead we compute the distance 
from the site to the curve. The walk must take at least this many steps before it can intersect 
the curve, so we can jump ahead this many steps in the walk before we check again for an 
intersection. 

A Proof of irreducibility 

In this appendix we prove that the pivot algorithm is irreducible in the half-plane and cut-plane 
that we have been considering. The proof is very similar to the proof for the full plane [H]. 
We show that for any self-avoiding walk in the half-plane (cut-plane, respectively), there is a 
sequence of pivot operations which "unfold" the walk into a straight line and such that each 
walk produced in this unfolding process is self-avoiding and remains in the half-plane (cut-plane, 
respectively). 

We first consider the half-plane. The restriction is that except for the starting point of the 
walk at the origin, the walk must remain strictly above the horizontal axis. We will show that 
the number of turns in the walk can be decreased by one. We denote the sites in the walk by 
a;(i) where i = 0,1, ■ ■ ■ , N. We will say there is a turn at u^i) if Lj{i — 1), a;(i) and a;(i -|- 1) are 
not CO- linear. 

We will consider cases based on the direction of the last step of the walk. If it is to the 
right, i.e., u{N) = uj{N — 1) + (1,0) we proceeds as follows. Let / be the largest integer such 
that the line y — x = I contains a site on the walk. So the walk is entirely below or on this line. 
Let i be the largest integer such that uj{i) is on this line. Since the last step of the walk is to 
the right, uj{N) is not on this line. So i < N. Since the first step of the walk in the half-plane 
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must be up, i cannot be 0. We take uj{i) as the pivot point and reflect the portion of the walk 
from uj{i) to uj{N) in the line y — x = I. The reflected portion of the walk lies entirely above 
the line, so the reflection does not produce self-intersections. Furthermore, since the walk was 
on or below the line, the reflection can only increase the y coordinates of points on the walk. 
So the new walk is still in the upper half plane. The walk before this reflection has a turn 
at uj{i) and the reflected walk does not. The reflection does not add any turns to the walk, 
so the total number of turns decreases by one. If the final step of the walk is to the left, i.e., 
uj{N) = uj{N — 1) — (1, 0), we use an analogous procedure with lines y + x = I to reduce the 
number of turns in the walk. 

Now suppose that the final step of the walk is either up or down, i.e., u^N) = u{N — 1) ± 
(0, 1). Consider the vertical line which contains this last step. First suppose that the walk lies 
entirely to the right of or on this vertical line. Let i < N he the largest integer such that there 
is a turn at uj{i). (Of course, if there are no turns the walk is a straight line and we are done.) 
The walk is a straight segment from uj{i) to uj{N) which lies on the vertical line. We take a;(i) 
as the pivot point and perform a rotation of 90 degrees (—90, respectively) if the last step of 
the walk is up (down, respectively). This rotates the segment from uj{i) to uj{N) to the left of 
the vertical line and removes the turn at uj{i). No new turns are added to the walk, so the total 
number of turns decreases by one. If the walk likes entirely to the right or on the vertical line 
containing the last step, an analogous argument shows the number of turns can be reduced by 
one. 

Now suppose that the walk contains sites on both sides of the vertical line which contains 
the last step of the walk. Let d be the horizontal width of the walk: 

d = max{x : {x,y) = uj^i), for some i,y} — min{x : {x,y) = uj{i), for some i,y} (38) 

We will show that d can be increased. Let I be the smallest integer such that the vertical line 
X = I contains sites in the walk. So the walk lies on or to the right of this line. Note that uj{N) 
is not on this line. Let i < iV be the largest integer such that u}{i) is on this line. We take oo{i) 
as the pivot point and reflect the walk from uj{i) to uj{N) in the line x = I. This increases the 
width of the walk. (The argument is the same as that given in fS[.) The reflection does not 
change the y-coordinate of points on the walk, so the new walk is still in the upper half-plane. 
Note that in the new walk the last step is in the same direction as before, i.e., either up or 
down. So we can repeat this procedure to increase d until we obtain a walk which lies entirely 
on or to one side of the vertical line containing the final step. When we reach such a walk we 
apply the procedure of the proceeding paragraph to reduce the number of turns by one. This 
completes the proof for the case of the half-plane. 

Now consider the cut-plane. The restriction now is that the walk cannot contain sites of 
the form (x, 0) with a; > 0, except for the starting point at the origin. We again consider cases 
based on the direction of the last step of the walk. If it is to the right, we proceed as in the 
half-plane algorithm. Note that the line involved, y — x = I, must have / > since the walk 
starts at the origin. (In fact / must be at least 1, but we do not need this.) The reflected 
portion of the walk will lie above this line while the cut, the non-negative real axis, lies below 
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it. So the reflection produces a walk that lies in the cut-plane. 

If the last step of the walk is to the left, a different procedure is needed to avoid producing 
a walk that intersects the cut. Consider the lines x — y = I and x + y = I. They intersect at 
(Z, 0) and divide the plane into four quadrants which we will describe as being left, right, above 
and below the point (Z, 0). We take I to be the smallest integer such that the sites on the walk 
lie in the quadrant to the left of (Z, 0) or on the lines. (/ is necessarily positive.) We then let i 
be the largest integer such that uj{i) is on one of the two lines. (It is not since the last step 
of the walk is to the left.) Note that uj{i) cannot be (/,0). We take uj{i) as the pivot point 
and reflect the walk from uj{i) to uj{N) in the line containing uj{i). The reflected portion of the 
walk will lie either in the quadrant above or below (/, 0) and so cannot intersect the cut. 

If the last step of the walk goes up or down we use the algorithm for the half-plane. There 
is a subtle point here. Recall than when the walk has points on both sides of the vertical line 
containing the last step, we chose / so that the walk is to the right of or on the vertical line 
X = I. For the half plane we could have chosen it so that "right" is replaced by "left." For the 
cut-plane this choice could result in a reflected walk that intersects the cut. To see that our 
choice does not produce a walk that intersects the cut, we observe that since the walk starts at 
the origin, / must be negative. The reflected portion of the walk will lie on or to the left of the 
line X = I, and so will not intersect the cut. This completes the proof for the cut-plane. 
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